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We theoretically study the polarization entanglement of photons generated by the biexciton cas- 
cade in a GaAs/InAs semiconductor quantum dot (QD), located in a nano cavity. A detailed analysis 
of the complex interplay between photon- and carrier coherences and phonons which occurs during 
the cascade allows us to clearly identify where the entanglement is generated and destroyed. A quan- 
tum state tomography is performed for varying exciton fine structure splittings. By constructing an 
effective multi-phonon Hamiltonian which couples the continuum of the wetting layer states to the 
QD we investigate the relaxation of the biexciton and exciton states. This consistently introduces a 
temperature dependence to the cascade. Considering typical Stranski-Karastanov grown QDs, for 
temperatures around 80 K the degree of entanglement starts to be affected by the dephasing of the 
exciton states and is ultimately lost above 120 K. 
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I. INTRODUCTION 

Coherent superpositions of quantum states allow to 
outperform the classical bit. Based on this feature the 
quantum bit (qubit) was introduced.^ In light of inte- 
grated quantum communication, natural candidates for 
its realization seem to be photonic systems. Not only 
do their two different polarizations serve as basis states 
for qubits in quantum computation, but also has encod- 
ing and manipulating of quantum information on photons 
made tremendous advancement: In recent years, photons 
as qubits have be en s ent successfully over fiber commu- 
nication channels) ^ quantum cryptography is already 
technically feasible^ and quantum teleportation of so- 
called entangled photon states was demonstrated Sup- 
ported is this development by the convenient fact that 
linear quantum optics is sufficient, when implementing 
quantum computing algorithms P 

Entanglement in its simplest form is a non-separable 
superposition of joint quantum states, in our case qubits, 
that show non- local quant um correlation. Among dif- 
ferent proposa kl 9 * 10 * 11 * 1 ^ ^ very promising solid state 
source for entangled photon pairs are semiconductor 
quantum dots (QDs) ™ In contrast to parametric down 
conversion,^ where entangled pairs are produced in a 
probabilistic manner with low efficiency (10~ 10 paramet- 
ric photons per pump photon 1 -) radiative decay of a biex- 
citon cascade in QDs provides an on-demand generation 
of entanglement, a crucial requirement for scalable quan- 
tum networks. 

The prospect of all-integrated photonic applications in 
combination with compact semiconductor devices raises 
the question of how robust and efficient an embedded 
entangled photon source is when subjected to losses and 
dephasing due to naturally occurring interaction with its 



surrounding host material. 

In this respect, Axt et al. investigated within a gen- 
crating functions approach in Ref. [T7] the dephasing of 
an exciton-biexciton-QD system, which is coupled to an 
arbitrary number of phonon modes and excited by a se- 
quence of 5-shaped classical light pulses, but did not con- 
sider the emission of single or entangled photons. With 
focus on photon entanglement Hohenester et al. consid- 
ered in Ref. [18] the elastic phonon scattering at the de- 
vice boundaries on a master equation level, assuming an 
asymmetry in the phonon coupling for different exciton 
states. 

In this paper we present microscopic calculations of a 
phonon-assisted biexciton cascade in an InAs QD embed- 
ded in a InGaAs wetting layer ( WL) ,™ The coupling of 
the discrete QD states to the WL continuum via multi- 
phonon processes^ leads to dephasing rates that signif- 
icantly limit the entanglement output efficiency for high 
temperatures (above 100 K). 

The coupled dynamics of occupation densities and 
photon-assisted states in the two-photon emission is 
treated with an equation of motion (EOM) approach, 
well-established in solid state optics to approach many- 
particle problems. This allows us to perform a quantum 
state tomography and thus calculate the concurrence of 
the polarization-entangled photons with varying, exper- 
imentally controlled external parameters like tempera- 
ture, exchange splitting, or WL carrier density. 

The paper is organized as follows: After the consid- 
ered system and calculated observables are introduced in 
more detail in Sec. [IT] the different interactions of the QD 
carriers with phonons and cavity photons are discussed 



in Sec. Ill Then the coupled equations of motion are de- 
rived, their dynamics solved and the results presented in 
Sec.[mT3]and[TV] 
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II. CHARACTERIZATION, GENERATION, 
AND MEASUREMENT OF ENTANGLED 
PHOTONS 

A. Entangled photons and qubits 

The two quantum states (|0) and |1)) that form a 
qubit \Q) can encode more information than a classi- 
cal bit, as a qubit generally exists in a superposition 
\Q) = ci |0) + C2 |1) with complex coefficients c\ and c^. 
In the model regarded here, the physical representation 
of the qubit states is the superposition of the horizontal 
\H) and vertical \V) polarization of photons. The true 
advantage of the quantum approach is reached when two 
photonic qubits Q = A, B interfere in a non-separable 
fashion and thus become ent angled P The wave function 
of any entangled qubit \tp}n can be written in the basis 
of the so-called Bell-statesP^ 
1 
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(\HV)±\VH)), 



where \HH) 



V2 
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= \H A H B ) = 
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\HV) , \VH) , \VV). Here, the entangled wave function 
cannot be expressed as a direct product \ip) q ^ \iPa)q ® 
\iPb)q of the wave function of each single qubit. Before 
going into details of our analysis, we note, that the final 
photon wave function generated by a biexciton cascade 
QD emission can be expressed in terms of only (see 
discussion below or Ref. IT4)l : 
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where is the phase difference between \HH) and |VV). 
Both systems are in a complete statistical mixture be- 
tween \H) and \V), but the coherence between \HH) 
and \VV) is kept. In the HH,HV,VH,VV basis, the 
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Following the Peres criterion J22l Eq. ([IJ indicates entan- 
glement since despite the remaining degree of freedom in 
6, the off-diagonal elements representing the coherence 
between the basis states can never be zero. If, however, 
these crucial matrix elements pyn — (W\ p pt \HH) arc 
zero, the system is said to have classical correlation only. 

Since we consider entangled photons generated by op- 
tical transitions in a biexciton cascade, we introduce the 
QD's electronic band structure and selection rules in the 
next subsection. 



B. Finestructure of quantum dots and generation 
of entangled photons 

Depending on the intensity, photo-excitation creates 
in direct gap semiconductors electrons in the WL con- 
duction band and holes in the WL valence band. The 
carriers subsequently relax into the QD and occupy the 
discrete energy shells following Pauli's principle Al- 
though relaxed, they still can interact with the WL via 
multi-phonon processes. To analyze entangled photon 
emission from a QD, its electronic eigenstates and se- 
lection rules for the light-matter interaction have to be 
known. 

Typically, when following a k • p approach, the band- 
structure of single-InAs QDs around the T-point can be 
approximately described by the two anti-binding s-states 
of the clcctronsand the six binding p-states of the holes P^l 
Due to spin-orbit coupling and strain effects the split- 
off (SO) and light-hole (LH) states lie energetically well 
below the heavy-hole (HH) state. Therefore, the rele- 
vant single-particle basis is constructed by the heavy- 
hole with total angular momentum = 3/2 and spin 
projection in growth-direction m^h = ±3/2 and the 
electron with J e = 1/2, m,j^ e — ±1/2, all shown in 
Fig. [TJa). The mutual Coulomb interaction will bind 
the carriers to electron-hole pairs and so lead to the for- 
mation of excitons. Depending on the configurations, 
given by the projections of the total angular momentum 
M = nij. e + rrij^h four exciton states arise, which can be 
characterized as optically inactive M — ±2 (parallel elec- 
tron and hole spin) and optically active M = ±1 (anti- 
parallel spins)) PS In a i?2d-symmetric QD these bright ex- 
citon states, denoted by \X±) are degenerated and couple 
to a± circular polarized light [here, cr + ( cr -) labels right- 
hand (left- han d) circular polarized light] as depicted in 
Fig.gb)™ ' 

The atom-like discrete energy levels of a QD can be 
pumped electrically or by photo-excitation.^ When the 
carriers occupy their lowest shell, two excitons can form a 
bound singlet state, a biexciton (\B)). Under emission of 
the first, <7_|_ or er_ polarized so-called biexciton-photon 
the system enters into an intermediate, optically allowed 
exciton state (|X±)). Subsequently, the QD relaxes into 
its ground state (|G)) by emitting a second, the exciton- 
photon. As a consequence of total angular momentum 
conservation both emitted photons are of opposite circu- 
lar polarizationj22l see Fig. fllb). Since the exciton states 
|X_|-) in symmetric D 2 d QDs are degenerate (i.e. no fine- 
structure splitting (FSS) 5 is present), the photons' de- 
cay path can only be determined by their polarization, 
otherwise they are indistinguishable. Thus, the cascade 
will produce a maximally polarization-entangled photon 
pair wave function |V , ) pt = (\HH) + \VV))/y2, which 
corresponds to the \$ + ) Bell-state.^ 

Although growth techniques of single QDs are very 
sophisticatedpS it is rarely possible to not have an asym- 
metry in the semiconductor crystal and so the non- 
classical correlation of the photons is often spoiled: Un- 
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(a) (b) D 2d /C 4v (c) C 2v 




Figure 1: (a) QD band structure with spin states s = |J;m 3 ). Relevant for the biexcitonic framework are only the electron and 
heavy hole (HH) states, (b) cascade of a symmetric QD with degenerate exciton states |X±). (c) cascade with a FSS 5 7^ 
(without biexcitonic shift). There are two possible paths. Either two photons with vertical V or horizontal H polarization are 
emitted into a cavity mode uiyju with FWHM k. 

For vanishing FSS S the which-path information is lost and the photons are polarization-entangled. With a FSS 5^0 the 
photons are not fully polarization-entangled. 



der strain the dots symmetry reduces to and the 
anisotropic electron-hole exchange interaction splits the 
exciton doublet into two states \Xh) = 1/ V%(\X+) + 
and \X V ) = l/V2(\X + )-\X-)i energetically sep- 
arated by the FSS S, shown in Fig. [ijc). These states 
couple to photon modes of orthogonal linear polarization 
along the direction of one crystallographic axis, labeled 
H and V respectively. This superimposes a which-path 
information onto the emitted photon frequencies and the 
degree of their entanglement is reduced. To efficiently 
collect the photons, the QD is placed inside a cavity sup- 
porting only two modes of different polarizations V, H 
with frequency ^ w^f v P^ These modes are assumed 
to be in resonance with the corresponding exciton-ground 
state transitions. Although energetically off-resonant, 
the biexciton-photon is emitted into the same mode, see 
Fig. [ljc). In principal, the exciton states can be tuned 
into near res onanc e again by applying an in-plane ex- 
ternal electricP^ 3 ^ or magnetic fieldj^ and indeed for a 
small FSS 8 generation of polarization-entangled photon 
pairs was demonstrated on a system operating at 10 K. 36 
Although the ideal case of zero splitting can be recov- 
ered, phonons as a decoherence mechanisms, in particular 
at elevated temperatures, will have an impact on the per- 
formance of a QD as a source of polarization-entangled 
photons. To provide a meaningful quantitative measure 
of the entanglement the next section will introduce the 
concurrence. 



C. Measure of entanglement — relevant quantities 

As shown in Eq. (JT|) , a measure for the degree of entan- 
glement is determined by the off-diagonal element in the 
polarization sub-space of the two-photon density matrix 
p pt , explicitly: 

Pvh := (VV\p pt \HH). (2) 

Quantum state tomography^ provides a measurement 
scheme which gives access to the different elements of 



p pt . They are experimentally reconstructed by measur- 
ing of the two-photon cross correlation function gfj(t, rf^ 
over a mean photon arrival time t. The function gfj(t, r) 
corresponds to the polarization correlation between a 
biexciton-photon emitted at time t and the subsequent 
exciton-photon at time t + rP^ 

4(t,r) <x (a\{t)a\{t + T)a 3 {t + T) a] {t)), (3) 

where i,j € {H,V}. The correlation function is written 
in terms of photon creation a\ and annihilation a opera- 
tors of the different photon modes cp. Fig.mc). We 
consider an experimental setup, where the time delay be- 
tween the two photons r is zero. This can be realiz ed by 
appropriately adjusting the distance to the detectorf 36 * 40 1 
The temporal dynamics of the corresponding density ma- 
trix element (ii\p pi \jj) is obtained when the second-order 
correlation function gfj(t, r = 0) is time averaged over the 
arrival times iP^ 

(ii\ P pt m(t) = \J*9U t '>°) dt '- ( 4 ) 

Thus, the source of entanglement can be rewritten as 
Phv oc ^ayaya H a H ^. A standard expression for the 
degree of entanglement is the concurrence 

C = 2\p HV \, (5) 

directly related to other accepted measures like the en- 
tanglement of formatiorP2 Ep or the tangle^ T (for ex- 
ample, C — VT). Here, C = 1 corresponds to maximum 
entanglement and C — to zero entanglement. To cal- 
culate the necessary dynamics of the expectation values, 

e.g. (^a\ / Oya H a H '^, a set of equations is derived in the 
next section. 
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III. MODELING AN EMBEDDED QUANTUM 
DOT AS A SEMICONDUCTOR SOURCE OF 
ENTANGLEMENT 

A. Hamiltonian and QD model 

The coupled dynamics of observables, such as Eq. |3]) 
can be generally derived from the system's Hamilton op- 
erator H via Hciscnbcrg's equation of motion: 



3D Bulk 



ihdt(d) = ([H,6] 



(6) 
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In this section, we discuss the used Hamiltonian to imple- 
ment our model system. Applying Eq. (|6j) to the photon 
correlation function gfj introduces a coupling to the elec- 
tronic degrees of freedom via the electron-photon inter- 
action. As discussed in Sec. |IIB| due to strong confine- 
ment, we assume that the QD single-particle states are 
energetically well separated and only the highest valance 
(v) HH state and lowest conduction (c) band s-shell form 
the biexcitonic framework, see again Fig. [T] In sec- 
ond quantization Fcrmionic operators describe the typ- 
ical electron-hole representation, where carriers are the 
heavy holes in the v band (operator h s ) and electrons in 
the c band (denoted by operator e s ). Here, the carrier 
spin state s = \J,rrij) is given for the holes (electrons) 
by T= |3/2;3/2) (|= |l/2;l/2» and |= |3/2; -3/2) 
(1= |l/2;-l/2». 

The complete Hamiltonian of a QD coupled to the WL 
continuum inside a nanocavity consists of various parts 
and is given as: 



H = H, 



QD,0 



TTC-C 



xVc-pt 

-"qd 



"QD.WL 



rrc-pn 
n WL,WL- 



(7) 



First, the kinetic energy of the confined QD carriers 
iiqp and their mutual carrier-carrier interaction Hq^ 
are introduced. The electron-hole pair in the QD does 
interact with the cavity photons of the quantized light 
field and so the free energy of the cavity photons 

Hq 1 is included as well. The free energy of semiconduc- 
tor bulk phonons Hq U and WL carriers appears, 
too. The interaction of the WL with the QD states via 
LO-phonons is considered in -ffq^wL anc ^ the electron- 
phonon coupling within the WL in ^wiTwl ■ The con t r i~ 
butions to the Hq part read: 



s 

Ho 1 =X^ L °^q' 



Figure 2: General system scheme. Two electron- hole pairs 
in a QD are coupled to the continuum of WL states via LO- 
phonon interaction pli,' 92 . The phonons are in a thermal bath. 
The QD has two levels with a conduction and a valence band. 
Here, the carriers interact with photons via the electron-light 
coupling elements M. The QD is assumed to be placed at a 
node position of the electromagnetic field in a cavity. Since 
its loss k = 10 [ieV is greater then the coupling strength to 
the field M = 1 ^teV, the system is in a weak coupling regime. 



The Bosonic longitudinal optical (LO) phonon creation 
(annihilation) operators at wave vector q are &1 (&_). 
Their dispersion is treated within the Einstein approxi- 
mation and Hlulo = 36 meV. Similar to the QD operators 



tu- 



rn 



ks 



) are creators (annihilators) of a hole carrier in 
the WL continuum of the valence band v with spin state 
s and wave vector k. For the WL carriers, we take into 
account only the hole contributions of the v band, moti- 
vated in the next section. The impact of spin-orbit cou- 
pling on the carrier's energy can be neglected in QDs 45 
and so, £j is assumed to be independent of the carrier's 
spin state. The general considered setup is displayed 
in Fig. [2] The next subsections discuss the remaining 
parts of the total Hamiltonian separately and in more de- 
tail. The electron-phonon interaction Hamiltonian leads 
to temperature dependent dephasing rates (see subsec- 
tion III All. The pure electronic part of non-interacting 



and Coulomb-interacting QD electrons Hq D + Hq^ is 
diagonalized and transformed into an excitonic basis. 
The new arising eigenvalues and eigenvectors will incor- 
porate the complete Coulomb contributions (see subsec- 
tion III A 2 1, which are energy shifts (e.g. ground state 
and biexciton shift) and the exchange splitting due to 
different spin combinations in the exciton states. In the 
excitonic basis, the electron-photon interaction is trans- 
formed and new orthogonal field modes arc derived (sec 
subsection IIII A 3 1 . 



Multi-phonon coupling of WL and QD 



H 



WL 



ks 

Off =^hwia\c 



Embedded in a host material, quantum confined elec- 
trons in Stranski-Krastanov grown InAs/GaAs QDs in- 
teract via LO-phonons with a continuum of two dimen- 
sional electronic WL states only a few ten meV away. 



5 



This leads to temperatu re dependent dephasing times for 
the QD state d 46 * 47 * 48 * 49 * and can be a first approach to the 
problem of temperature dependent generation of entan- 
gled photon pairs. On the time scale of several nanosec- 
onds regarded here, longitudinal acoustical phonon^ 50 * 51 * 
only have a minor impact and are not considered. 

Depending on the dot sizej^ the QD valence band is 
typically more than one, but less than two LO-phonon 
energies fkJho energetically separated from the WL band 
edge, see Fig. [3] To effectively connect the QD states 
with the WL, up to two-phonon processes have to be 
taken into account. Within the two-phonon limit the 
influence of the WL conduction band on the QD elec- 
trons can be neglected, because here, more than two LO- 
phonons are necessary to bridge the energy gap to the 
WL and therefore the dephasing is determined by the 
hole-WL interaction. Moreover, the Coulomb interaction 
of the WL carriers is not included as the carrier densi- 
ties considered here are lowP^ Under these assumptions, 
microscopic dephasing rates are derived by using an ef- 
fective Hamiltonian approach, which originates from a 
multi-photon theory.^ From the Hamiltonian in Eq. ([T]) 
the following parts contribute to the LO-phonon induced 
dephasing: 



H. 



(>D,WL — ^QD,0 

+i? pn + 



rVc-pn 
n QD,WL " 



rrc-pn 
fl WL,WL- 



(8) 



The phonon mediated interaction between QD holes and 
WL states and the carrier-phonon interaction within the 
WL are given by 



rrc-pn 
■"QD,WL 



= EE 90k.H™ks(h + *>-q) + h ' a -: 



kq 



#wl"wl = EE 9kw^lsK's(K + fr-q) + h.a.. 

s kk'q 

The Frohlich coupling elements are g q k , which can be 
found in Ref. [MJ Within a projection operator based 
theoryMl Eq. (jgj) j s mapped onto the resonant WL states 
only and becomes: 



QD.WL 



= H, 



QD : 



WL,0 



H, 



ett: 



(9) 



with the effective LO-phonon-assisted WL influence on 
the QD holes in H ff. All other off-resonant contributions 
are implicitly included in the coupling elements of the 
effective Hamiltonian. Taking only two-phonon processes 
into account, H c g reads 



E E 

s q 1 q 2 k r , 



0„}l, 3 h a 1 



Pwc,s W \r r 



+E E 

s qi q 2 k res 

with the effective coupling elements 

2-^i -WL ^QD %l , 

k £ vk ~ £ y - n ^LO 

k^k ros 



r,qiq2 



(10) 
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Figure 3: Wetting layer system scheme. Seen in blue is the 
2D dispersion of the WL holes ft 2 k 2 /(2ml). Marked on the 
energy scale are (i) the WL band edge £^ L which is ener- 
getically separated from the QD v shell e? D by more then 



a single-phonon energy: Ae^ D 



-QD 



> hojLO, (ii) 



the WL states at k ros resonant with a two-phonon transition 
e3£ M = e?° + 2^lo to the QD. 

The transition from the resonant states to the QD pass 
through intermediate states at eJX L with probability ampli- 
tude Sk 2 k, • 



They contain Pauli blocking terms and therefore de- 
pend on temperature and WL carrier occupation (1 — 

(^ks^ks))-^ Th° WL holes in H e g have an energy ex- 
actly two phonon energies away from the QD state en- 
ergy. A transition from these resonant WL holes at e^. L 
to the QD shell takes place under simultaneous emission 
of two phonons. The whole process is energy conserving: 



-WL 
'fk rc 



£ QD _ 2/kjLO- Within time-energy uncertainty 
carriers relax by a higher-order Markov process. Here, 
in the transition to the intermediate state at e^ L energy 
conservation is violated, since the hole state at k is less 
than Swlo from and more than fio>LO from the QD 

state energetically separated, see Fig. [3] The probability 
amplitude for the intermediate transitions are g q k and 

.<aV 



Equation (10 1 shows that all possible transitions 



between QD and WL are mediated by all off-resonant 
WL states k. The strength of the coupling element is 
determined by to what extent the energy conservation 
condition in the denominator is met in every phonon- 
assisted electronic transition. We can use the effective 
Hamiltonian Eq. Q to derive relaxation and dephasing 
rates using Heisenberg equations of motion, where the hi- 
erarchy problem is treated within a born factorization.^ 
The calculations lead to the following equations for the 
QD states: 



dt 
d_ 

dt 



(e\h\) = -(ihn + r w , s )(ethi), (ii) 

(hth s ) = -2T w , s {htk), (12) 



with the QD gap energy KIq — — e^ D and the WL- 
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Figure 4: Plot of the phonon-induced relaxation rate I\„ as a 
function of temperature T. At about 80 K T w has approached 
1 ns _1 and starts to contribute strongly to the decay of the 
QD states. 

induced damping rateP^ 

T w , s = [(1 - / M )(n + l) 2 + fh, s n 2 ] j s . 
The damping 7 S is given by: 



(13) 



7s = 



J f d 3 9l d 3 q 2 (pS# + Plti ) ■ 



(14) 



In Eq. ([13]) f h , s = Sk rca (^k rcBS ^k rcsS ) is used for tne WL 
hole density at the resonant energy, which in the carrier 
low-density limit is assumed to be zero. Note, that this 
implies T w — T w ^ = F w ,l- Since the system has relaxed 
into a quasi-equilibrium, the phonon bath is described 
by the Bose-Einstein distribution n = (b^b). Figure Q 
displays the temperature dependence of T^, . In the next 
section an excitonic basis will be introduced. The damp- 
ing rates in Eqs. ( 11 ) and ( 12 1 will be used in this basis 



to account for the relaxation of the QD electrons.^ 



2. Carrier-carrier interaction and exciton representation 

The Hamiltonian in Eq. ^ accounting for the QD 
carriers and their interaction via the Coulomb potential 
-Hcoui = ^qd o + QD i s conveniently rewritten asP^ 



hlk 



(15) 



V™e\e s e\,e sl+ V hh h\hji\;hs- 

ss' 

+ 2V?f, efh\h s ,e s , - 2V he e&&+,A.,) . 



Table I: Numerical parameters. 



Parameter 


Symbol 


Value 


electron effective mass 


m e 


0.043771^ 


hole effective mass 


m h 


r\ a rr» 16 1 1 


LO-phonon energy 


fru>LO 


r>^ a T 46 1 1 

36.4 meV^ 


QD band gap 




1.5 eV 


hole binding energy 




1 Fi flCjJi n 


Coulomb parameters 




25 /ieV 




t rex _ i/ex 


fieV 


photon lifetime in a cavity 


K 


10 AteV 


electron-photon coupling 


M 


1 ^eV 


photon lifetime 


Trad 


50 ps" 1 



The first term contains the non-interacting electrons with 
gap energy hflo- The second term accounts for the re- 
pulsion of carriers within the same band, whereas the 
last term gives attractive direct Coulomb interaction 
and repulsive exchange interaction between carriers in 
different bands. The corresponding Coulomb elements 
V ee , V hh , V he , mediate the interaction. Responsible 
for the fine structure splitting 6, compare Fig. [l] is the 
exchange splitting V*?, which describes the repulsion and 
attraction forces induced by different spin-conformations 
of electrons and holes. As shown in App. [A] the FSS can 
be expressed by 



8 = 2\Vfi\. (16) 

To simplify the notation we will refer to Vf* as V cx . 

In principle, all matrix elements of the Coulomb in- 
teraction can be microscopically calcul ated, whe n the 
single-particle wave functions are known E2ESE1E2 How- 
ever, their values only have a quantitative impact on the 
results. Therefore, within a reasonable range, they are 
used as model parameters, measured in experiments 
Using the space spanned by the new exciton operators, 
derived in App. [a] the Hamiltonian is rewritten asP^ 



H Coul = huj G G^G + kujH X' H X H 

■>B 



huj v Xlx v + hujR B^B 



with G the ground state, X H / V the exciton and B the 
biexciton annihilation operator, corresponding to the ex- 
citonic level structure as depicted in Fig. [TJc). Note, 
that within this diagonal representation the derivation of 
equations of motion via Eq. © is trivial for the Coulomb 
interaction, as operators of different states commute. 



3. QD electron-photon interaction 

Commonly, the Hamiltonian of the electron-photon in- 
teraction is taken in rotating-wave approximation 



h.a., (17) 
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Figure 5: Equations of motion, truncated to the pure cascade 
scheme. The dark blue quantities represent densities which 
do not contribute to entanglement, whereas the dark orange 
and red quantities directly generate a crossing of the different 
paths in the light orange boxes and are crucial to entangle- 
ment. 



where the corresponding spin states couple to circular 
<t± polarized light (see Sec. IIB| for more detail), i is the 
mode of the emitted photon and M denotes the electron- 
photon coupling matrix elements. The Hamiltonian of 
the light-electron interaction is expressed with the exci- 
ton operators, derived in the App. 15] 



= hM (&X H a\ H + &X V aU (18) 



a] v ) +h.a.. 



We assume that our QD is placed at a node position in- 
side a nano-cavity, which provides two different modes for 
the different polarizations H and V corresponding to dif- 
ferent frequencies ujh 7^ wy Since only two modes exist 
within the cavity we investigate a cavity-enhanced biex- 
citon cascade. However, we remain in the weak coupling 
regime since the cavity loss n = 10 /ieV is greater then 
the coupling strength to the field M — 1 /ieV. Regard- 
ing only these two modes, the electron- light interaction 
Hamiltonian can be written in a compact form 



h$ = m 



H H 



-X%B tf H 



v u v 



X y B a v 



h.a. 



(19) 



At this status, the total Hamiltonian is written with the 
new exciton and photon operators (iJ, V) . To discuss the 
polarization entanglement between two emitted photons 
we proceed and determine equations of motion, which 
govern the relevant observables' dynamics. 



B. Equations of motion 

To derive the coupled dynamics of photon-polarization 
coherences and electronic transitions with Heisenberg's 



equation of motion Eq. ([6|,l 65 l 66 lthe total excitonic Hamil- 
tonian 

H = J2,- H v ^fStfij + huj G &G + hw H x\jX H 
+ huo v X V X V + huo B &B + hM (&^X„a}„ 
+ &X v a v + X\jM H - X v Ba v + h.. 



..a. 



in conjunction with the temperature dependent relax- 
ation rates T w given in Eq. ( 13 1 is used. The latter lead to 



an additional decay of the QD populations and dephas- 
ing contributions to the QD transitions. Those contribu- 
tions are derived via the effective Hamiltonian approach 
in Sec. |III A 1| and are consistently included by higher- 
order Markov approximations of the phonon-mediated 
interaction between carriers in the WL and QD.^SJ 

An overview of the, at a first glance complicated cou- 
pled dynamics of the considered correlation functions 
(derived by Eq. [fjj, is given in Fig. [5] Going through the 
scheme, step by step we unravel the consequential inter- 
play of the different quantities. An initially given biexci- 
ton density can decay via two possible paths (left B, right 
V) and generate a photon pair Pij(i,j € {H, V}). In a 
first step, a photon-assisted coherence builds up x\B a\ 
(light blue box), which then contributes to (i) a cross- 



polarization coherence X V X H a v a H (red box) partic- 
ularly important to achieve entanglement in pvh- (h) 
a two-photon coherence & B a\a\ (light orange box), 
which also leads to an interference of the two path and 
thus contributes to pvh- (hi) a combined exciton-photon 
density X\X i aj^ (dark blue box), which does not influ- 
ence the degree of entanglement. This gives meaningful 
insights to the underlying physics and origin of polariza- 
tion entanglement. Note, that we are in a weak coupling 
regime and only spontaneous emission in the cascade is 
taken into account. 

The concurrence C as a measure for the degree of entan- 
glement is determined by the photon density matrix, cf. 
Eq. (I5]|, and defined via the off-diagonal element pvh '■ 



d t (a v a v a H a H > = 



(20) 



2i (u^ av - lu c £ v 



- 2ire) (^a v aya 1 :(i n 



+2iM U&X H a v a 



H 



Beside its damping due to cavity losses^ chosen to be 
k = 10 /UeV, the two-photon correlation pvn is driven 
by two higher-order quantities. Both include an exciton- 
ground state transition under emission of a photon of 
opposite polarization as the state would allow, e.g. 
G^X H a v (the complex conjugate of &X v a H ). The 
transition process takes place under presence of a pho- 
ton coherence a v a H generated by the previous biexciton- 
exciton decay, see light-orange boxes in Fig. [5] As these 
terms already include a single-photon coherence and gen- 
erate a second one leading to a two-photon coherence, 
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Figure 6: Temporal evolution of the non-integrated off- 
diagonal two-photon density matrix elements at T — K 
(real (solid) and imaginary (dashed) part) for V ex = 1 peV 
(orange) and V ex = 10 fieV (red). With increasing V ex the 
oscillations become more rapid. 



they are exactly the terms one would expect to contribute 



to pvh- For a small FSS, the fixed cavity frequencies u> 



and 



cav 
V 



are in near-resonance and pvh will slowly os- 



cillate on the timescale given by the corresponding FSS, 
see the orange curves (all at T = K) in Fig. [6] For 
an increasing FSS on the other hand both frequencies 
are detuned and pvh shows a strong oscillating behav- 
ior, compare red curves in Fig. [6] Here, the temporal 
mean value of pvh is close to zero and thus no entan- 
glement in a measurement is observed, compare with red 
curve (all at T — K) in Fig. [7] for the integrated pvn- 
In a physical interpretation that means the two different 
decay paths are distinguishable, so either the photons 
are entirely emitted in the H or V cascade, but there 
is no overlap which is only generated by transitions like 



G^X 



containing both V, H indices. The which-path 



information is conserved. If there is an uncertainty in the 
decay path, the photons become partially polarization- 
entangled. 

The quantities in this section are damped by the LO- 
phonon-assisted WL influence as motivated in section 
The calculated damping rates T w from Eqs. ( 11 ) 



IIIA1 



and ( 12 1 correspond to a Ti-time and are incorporated 



like the radiative dephasing r ra( j in the Weisskopf-Wigner 
theoryPSES Both occur and lead to an overall damping 
of r = r rad + T w : 



d t {B^BJ (x -4T {B^BJ , 
dt^xjBa^ST^XjBa, 

The exciton Xi is damped by the mere presence of the 








500 1000 
Time [ps] 



1500 



Figure 7: Integrated two-photon matrix elements at T = K. 
Their steady-state values give the quantum state tomography. 



empty WL states as they disturb the system and lead 
to a decoherence. This introduces a temperature depen- 
dence to the cascade, which is inherit to the perturbation 
induced by the coupling to the WL and thus included in 

r„ 



T W (T) motivated in Sec. IIIA1 



By numer ically solving the set of Eqs. (|20j p]|cT4| 



given in App. [Cjwe will investigate how the temperature 
affects the concurrence. Later, the complete temporal 
dynamics of the cascade presented already partially in 
this section is resolved. 



IV. RESULTS - DYNAMICS OF THE 
BIEXCITON CASCADE AND QUALITY OF 
ENTANGLEMENT 

Since p pt is experimentally reconstructed by photo- 
counting experiments all elements of p pt are given by 
time-averagingPS Recalling and employing Eq. ([3| to the 
results of Fig. [6] where already the temporal evolution 
of p pt elements is given we get its integrated elements. 
As can be seen in Fig. [7J the diagonal elements have a 
continuous positive slope and start to saturate around 
0.5 ns. The steady state which determines the quantum 
state tomography is reached at 1 ns. However, the sit- 
uation is very different for the off-diagonal elements, as 
they are complex quantities that oscillate when V cx =/= 0. 
Its absolute value (important for the concurrence) shows 
a non-monotonous behavior in Fig. |jj Obviously, the 
concurrence is lost for a FSS higher than V ex — 10 p,eV. 

Since all steady-state elements of p pt are given by time- 
averaging, C will drop for increasing V cx as the inte- 
grated pvn did. This can be clearly seen in the quantum 
state tomography shown in Fig. [8] The diagonal and 
off-diagonal contributions are still in the same order of 
magnitude for V cx — 1 /ieV cp. Fig. |8ja), but a loss can 
already be seen. For a larger splitting pvh vanishes in 
Fig. [8|L). Figure [9] constitutes the central result of this 
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(a) Vex = 1 A*eV 
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(b) V ex = lO^ieV 
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Figure 8: Quantum state tomography for (a) V G x = 1 AteV 
and (b) V cx = 10 /ieV. 
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Figure 9: Plot of the concurrence, (a) displays C at K 
as a function of V cx only, (b) illustrates the influence of the 
temperature due to phonon-induced WL influence at Kx = 
/xeV. Both temperature and V cx dependence are shown in 
(c). 



work - a temperature dependent study of entanglement. 
First, Fig. ]9ja) shows how the entanglement is lost with 
increasing FFS as a continuous function of V ex - Here, 
the FWHM is determined by the values of the Coulomb 
parameters. When temperature effects of the WL states 
are taken into account, the concurrence can be spoiled 
even in the ideal situation of degenerate exciton states, 
see Fig. |9jb) . For low temperatures C will remain unaf- 
fected by the WL-induced dephasing T, since the scatter- 
ing times are well above 1 ns, cp. Fig. [4] Starting at ap- 
proximately 80 K the WL starts to affect C as T reaches 
1 ns -1 , which corresponds to an energy of 0.7 /j,eV close 
to the optical coupling strength of M = 1 fieV. The en- 
tanglement decreases for zero V cx until it is entirely lost 
for temperatures beyond 100 K. For a higher FSS with 
V^ x 7^ 0, Fig.^c) shows, that the degree of entanglement 
is lost slightly earlier around 80 K. 

Finally, to pick up on the topic of temporal dynamics 
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Figure 10: Dynamics of the cascade biexciton — > exciton- 
and-one-photon — > ground-state-and-two-photons. Results 

for vanishing V cx = fieV at T = K. (^X^ Xho^jj^H^ an d 

Phh are enlarged by a factor of 10 4 and 10 5 , respectively. 
The inset shows the dynamics on a logarithmic scale. 



of the cascade already addressed in Sec. |III B] let us con- 
sider a direct, single path leading to no entanglement. 
Therefore, we will follow the blue HH (left) path in 
Fig. [5] The biexciton density (B^B) decays exponentially 
with 4 r giving rise to an intermediate coupled exciton- 

photon state (&HX H a\ja H \. Subsequently, when this 
state is sufficiently populated it decays under emission of 
the exciton-photon and a two-photon density phh builds 
up. In the given range of parameters (see table]!]), Fig. 10 
shows that the decay cascade happens on a ps time scale. 
Even at low temperatures and V cx = fJ-eV, due to high 
cavity loss n (compared to the optical coupling strength 

M) both (x\jX H a} H a H \ and phh arc only weakly oc- 
cupied. The inset in Fig. 10 is a logarithmic plot of the 
dynamics which clearly shows the different lifetimes of 
the involved quantities. 



V. CONCLUSION 

In summary, we showed that - based on a Heisenberg- 
equation approach - the density matrix of the photon- 
polarization subspace of a biexciton, all intermediate oc- 
curring states of the cascade as well as their dynamics 
can be microscopically calculated. The interaction of the 
dot states with the WL via LO-phonons was included 
within this approach and gave rise to a strong reduction 
of the concurrence for temperatures above 100 K for typ- 
ical InGaAs self-assembled QDs. 

The diagonal interaction of the QD states with lon- 
gitudinal phonons is a major contribution to dephasing 
effects psESI which will ultimately influence the quality of 
entanglement 
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Our conclusion is, that regardless of their impact, the 
inherit coupling to the WL imposes another fundamen- 
tal limit to high temperature generation of polarisation- 
entangled photons in solid state devices. 



which are not eigenstates of f/coui- For convenience, via 
solving a diagonalization problem, new exciton operators 
X H / V are introduced. Here, H and V refer to the lin- 
ear polarization of the emitted photons in the biexciton 
cascade: 
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Appendix A: CARRIER-CARRIER 
INTERACTION AND EXCITON 
REPRESENTATION 

When constructing the single-particle basis which is 
to be diagonalized we can employ the fact that only a 
fraction of all possible QD states will contribute. First, 
not all transitions are optically active. A transition of 
a conduction band electron under emission of a photon 
must conserve the angular momentum. The electron- 
photon matrix element leads to selection rules that have 
to be obeyed and so only electron-hole pairs with opposite 
spins couple. Second, and more important, the Pauli- 
Principle forbids two carriers to be in the same state 
(that is the spin state s = \J,irij)). 

With these considerations, only four states remain to 
determine the system dynamics, defined by the following 
operators: 



G 



/ij^ /i^ /i-j- 



Ground state operator, (Al) 
Biexciton state operator, (A2) 
Exciton state operator 1, (A3) 
Exciton state operator 2. (A4) 

Commuting these operators with Hcoui = -^qd o + ^c-o 
we note, that ground state G and biexciton operator B 
are already diagonal: 



HuiG G, 



H, 



Coul 



B 



huj B B, (A5) 



with Hlog = 0, choosing zero as ground state energy, 
and hu B = 4E C + 2hfl + + Vff with E c = 
I (y cc _|_ y vv — 2V VC ). Investigating the exciton opera- 
tors X±, we find: 



H, 



Coul, 



, x+ = (e c + nn Q + v?f)x+ - v?f x_ (A6) 



Hcoui,*- = (E c + tXlo + Vff)X- -Vff X+. (A7) 



The exciton operators in Eqs. (A3) and (A4| emit into 
circular polarized light modes (er + and ct_), if the non- 
diagonal element is zero. In a reduced C^v symme- 
try of strained dots, the off-diagonal element V** is non- 
zero and leads to a superposition of the exciton states, 



X H =u^X + + v% X_ , X v =u\X + + u\X_ . (A8) 
The unitary transformation coefficients are given by 



,H 



-(l + A*)-*A e , 



(1 + A^)- 



where A e = V m (Ec + t&l +Vft -T^h)" 1 = -V CX (E C + 
HQ + Vff -Hujv)- 1 and V cx = Vff = Vff. Within the 
new excitonic basis the electron operators (G^\ Xjj , 
Xy \ B^) are eigenvectors to #qd.o+^qd- Their corre- 
sponding eigenvalues are (hujQ,hw}j 1 hioy^ ! hujB) with ex- 
citon energies Hloh and hujy in the two-particle basis: 



ii v = 2 (- / '<" 



2hfi Q + V\ 



n 



Fff ± 5) 



(A9) 



Responsible for the fine structure splitting 6, compare 
Fig.[T| is the exchange splitting V™, which describes the 
repulsion and attraction forces induced by different spin- 
conformations of electrons and holes. In the most general 
case, the exciton states could differ in energy due to con- 



tributions like V. 



TT 



and VjJ. 



Given these elements, the 



most general fine structure splitting can now be expressed 
quantitatively as S := 



T/ex _ T/cx 

TT U 



■4|^|2. 



How- 



ever, it is reasonable to assume that in semiconductor 
QD no spin-preferences exist, thus V££ 



leads to 



Vff = 0, which 



s = 2V n\ 



(A10) 



H 



In our case of no spin-preferences, where {Vff—Vff) = 
it follows that A e = —1, u 
and thus explicitly 

1 



l 

V2> 



x, 



x_ 



XV ~T2 



x+ - x_ 



(All) 



Appendix B: QD ELECTRON-PHOTON 
INTERACTION 



Starting with Eq. (17 1, we now switch to the new exci- 



ton operators X H iy by inserting the unity relation of the 
electron-hole picture into Eq. (171. After normal ordering 



and using the two-electron assumption, the electron- light 
interaction can be written asP^ 



H$ = HM 



X-B\a 



it 



(Bl) 



hMj2 (&X_ + X+B\ a\ a _ + h 



11 



The electron-light interaction Hamiltonian -Hqd* is trans- 
formed, when the exciton operators are replaced with: 

X + =-L (x H + X V ) , X- =i - *v) ■ (B2) 

It is convenient to define new photon operators: 



1 / 



t J_ At 



1 A 



1 At 
+ 



(B3) 



The Hamiltonian now takes the form: 
^QD = KM J2 (&X H a\ H + &X V a\ v ) (B4) 



dL- + h.a. 



Appendix C: EQUATIONS OF MOTION 



The temporal evolution of the driving terms in Eq. ( 20 ) 
is given by 

d t (G ] X H a} v a} v a H ^ = 

i ( — u>h + 2c^ av - u c h v + ir + 3i/c) (&X H a f v a ! v d,j^ 
-2iM ^4i H fl{,a ff J + iM /(^.Ba^a^ (CI) 

and 

d t (XyG a v a H a H J = 
i (lo v + cu^ v ~ 2uf£ v + ir + 3i«) I XyG a v a H a H ' 



2iM ( X v X H a v a H 



iM(B t G 



(C2) 



The driving terms of the two-photon density matrix in 
turn couple to combined exciton- and photon coherences 
X v X H a v a H and to the direct decay channel from \B) 
to \G) emitting two photons with the same polarization 
& BcLyOy., see orange box in Fig. [5] Crucial for en- 
tangling the two decay paths is the exciton coherence, 
assisted by a photon coherence, see red box in Fig. [5] 

d t (x v X H a v a H ) = 



i (lo v ~lu h + cjy* v - + 2ir + 2i«) (x v X H a v a H 
iM {G^X H a v aya H ^ +\M (B^X H a H \ 

(C3) 



-iM ( X V G a v a H a H ) + iM (X v Ba v 



't OAt 



In this equation, the two paths interfere. The influence in 
the two-particle correlation (x v X H a v a H \ increases the 



degree of entanglement as this term couples back to the 
driving terms of pvh, Eq. (CI ) and ( C2 1 . Here again the 



resonance condition of the frequencies is essential (ujy — 
J W = A high detuning 8 will diminish 



the contribution of Eq. ( C3 1 to the cascade and both 
paths cannot interfere. 

The other characteristic and important quantity in the 
two-electron biexciton-cascade situation (cp. with two 
coupled are the two-photon polarizations 



dt I &Ba v a v ) = 



i(-w fl + 2w^ av + 2ir + 2i«) (&Ba 



l V V 



(C4) 



+iM (G^X H a v a v a H \ -iM Xya v a v ay) 



"2 LI/ (XyBOy 



and 



d t (B^G a H a H ) = 



i{w B - 2w# v + 2ir + 2i«) (B^G a H a H ) (C5) 



AM (X f H G a} H a H a H ) +iM lXyGa v a H a H 



2iM (B^X H a H 



(C6) 



Each path in the cascade has one biexciton-to-ground 
state transition like & Ba v ay. Its dynamics cou- 
ples the biexciton-to-exciton transition X v Ba v with 
both exciton-to-ground state transitions G'Xi. Remark- 
ably, the origin of the entanglement is directly visible, 
since a quantity of a different path enters in Eq. (C4|: 

(G^X H a v a v a H 



Here again, the two paths interfere. 
For maximum entanglement the contributions of the dif- 
ferent paths & Xh and G'X V to the expectation val- 
ues should be equally weighted. The photon-assisted 
biexciton-to-exciton transition enters in the two-photon 
polarization and drives this quantity via the biexciton 
decay: 



dti&XvV 



H H 



i (uj b -UJH- u c h V + 3ir + i«) {B^X H a H 
iM (x\jX H ($ H a H )+iM (x v X H a v a H 
\m(b*B)+\M (&G a H a H \ , 



(C7) 



dt(X v Ba v ) = 
i (-u B +0Jv + 4 av + 3ir + in) (x^My 



(C8) 



iM (B^B 



-iM (&Ba v a v 



-iM (XlX H ala H )-iM (xlx v a\,a v ) 
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The occurring biexciton as well as the intermediate 
exciton-photon densities are driven by the biexciton- 

exciton transition ( XjBa, 



(C9) 



d t ^ H X H a) H a H 
- (2r + 2k) (xj 1 X H a' H a H 



-2 Tin ! .1/ (xjjBa^ 



M ( XjjG alj a H a H 



To complete the set of equation, two higher-order photon- 
assisted exciton-to-ground state transitions of the direct 
and thus not entangled path are necessary: 



d t (&X H a 1 H a i H a H ^ = 

i (-loh + lu c h v + ir + 3i«) (d^X^a^a^ (C12) 
-2iM /x i H X H a i s a H '\+iM ^Ba^a^, 



- (2r + 2k) (x v X v a v a v 



(CIO) 



-2 Till ( 1/ [XyBOy 



— M (XyG a\ra v a v 



From the perspective of the cascade our course of action 
so far put the cart before the horse since the actual dy- 
namics start with a loaded biexciton density (b^ B^j. In 

the visualization of the complex interplay, Fig. [5j we fol- 
lowed a bottom-to-top trail through the cascade, starting 
with the concurrence determining pyn- The biexciton 
(^B^ B^J as the top element of the scheme decays via the 
H or the V intermediate exciton-to-ground-state path 



dt(B^B) = -4T(BTB 



(Gil) 



-2 Tm ! .1/ ( X^Ba^ ) - M (x v Ba v 



d t (d ] X v alala v ^ = 

i ( — U)y + + IF + 3i«) (G^ XyOyOyhy^ 



2iM (XyXyOyhy 



iM (&Ba\,a\ 



(C13) 



With these polarization Eq. ( C12||C13 l, the diagonal el- 
ements i = H, V of the density matrix of the polarization 
subspace are given, too: 



-4k (alala^t 



— 4 Im ( M I &X]a]a 



(C14) 
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